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SUMMARY 


omaaM? ' 


The errors involved in using the perfect-gas lav pv = RT the 

assumption of constant heat capacities are evaluated. The Basic equations 
of gas flows taEing into account these phenomena separately htu^ at the 
same time are presented. 


INTRODUCTiaN 


The conventional method of ohtainlng higji Mach nuniberB for 
aerodynamic tests is to accelerate the air by means of a pressure 
difference so that the random kinetic energy of the molecules of air at 
rest is converted into kinetic energy in the test section. For very 
hi^ Mach numbers this may occasion high stagnation temperatures nrifl 
pressures which introduce effects due to the vibrational heat capacity 
and molecular forces and size such that the perfect-gas law pv = ET 
and the assumption of constant hean capacities may be no longer 
sufficiently accurate to evaluate gas flows. 

It is the pjurpose of this paper to present formulas which are 
suitable for handling such problems and to point out the magnitude of 
the errors that may be involved in using the perfect-gas law snd the 
assumption of constant heat capacities. 

Tsien (reference l) has published a theoretical discussion of 
this problem in which certain approximations were introduced in order 
to obtain solutions that were in a very neat form when the Irperfect— 
gas effects were moderate. A ccanparlson of Tsien’ s results with 
this work is presented to show the magnitude of these approximations. 

In England Goldstein has previously investigated this problem 
at moderate temperatures and pressures in order to prove the small 
magnitude of imperfect-gas and vlbratlonEil— heat— capacity effects in 
most supersonic wind tunnels. This report indicates, in general, the 
range in which these effects are small but' does not present formulas for 
hand ling problems in gas dynamics when these effects are large. 
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Mie present paper Is arranged in tiie following tiiree parts: 
temperature effects on perfect-gas flows due to variation of Weat 
capacities; in^jerfeot-gas effects on gases without variation of 
heat_ capacities; and. gas -flows in which, hoth effects are present. 

This is done since idle formulas in the first part may prove use- 
ful to those dealing with the flow of hot exhaust gases a-Tui since 
it may bring out more clearly the differences between the two effects. 


SYMBOIS 

a term in Yan der Waals* equation correcting for the 

effect of molecular farces 

b term in Yan der Waals* eq,u!itian correcting for the 

effect of molecular size 

c speed of sound, feet per second 

Cp heat capacity at constant ]>ressure 

Cy heat capacity at constant volume 

E , energy, foot-pounds 

M , Mach number (w/c) 

p pressure,' pounds per square foot 

E gas conetant 

T absolute teirrperature, degr«>es Eahrenheit 

V specific volume, cubic feet per slug 

w velocity, feet per second 

7 ratio of heat capacities (C^/C^) 

p density, slugs per cubic foot 

9 characteristic temperatture of molecular vibration 

Subscripts: 

o B'bagnation coEBiitions 


c 


criticeil conditions 



NAGA EM Wo. l8Jl4 


3 


Errars Involved in Aseuming Constant Specific Heats in the 

Presence of High. TemB«ratures in a Perfect Gas 

For a perfect gas with constant heat capacities the eguatlon for 
conservation of energy of a steady isentropic process may he "written as 


CpT + |»r2 = 


If this 6 q,uation is combined vlth the equation for the isentropic speed 
of sound c 2 a 7 ET the resulting equation is 



( 1 ) 


If the eipanslon is Isentropic, the pressure and density ratios 
corresponding to the Ifech number are 


and 



( 2 ) 


However, if the temperature of the gas is hl^ enough the heat 
capacities may not he assumed constant because the vihrational degrees 
of freedom of polyatomic molecules are ezcited. The variation of the 
equllihrium value of the heat capacity at constant voliuns of a perfect 
diatomic gas is found from qtiantum mechanical considerations to he of 
the form 


E 2 \tJ 


( 3 ) 


■where 9 is a constant depending on the gas. The formula may he used 
for the mlsrbure air if the value of 9 is placed equal to 5526 when 
absolute temperature is measured in degrees Fahrenheit. (See the 
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appendix.) The value of- the heat capacity at; constant ^pressure for a 
perfect gas Is then 


^ = I + fif 

E"2 W te/r_3.)2 


(h) 


Figure 1 Is a plot of equation (h) aod shows that the heat capacity 
may not he considered constant ahove 600 ° F absolute. 

When the heat capacity at constant piressure varies according to 
equation (i)-) the energy equation must he ‘irrltten 


E 



^ dT H- 


0 


(5) 


Substituting equation 


(!<■) Into equation ( 5 ) and Integrating yields 


e®/^ - 1 


+ gR0 

_ 1 


(S) 


The Mach number Is obtained from equation (6) by dividing through 
by , 7 ET = c^, which gives 


where 




®A'o _ 1 e®/T — 1/ 


Cp/R _ 7(e®/^ - if + 2(§)V/^ 
(V/R 5 (^0/T _ + 2@ )^e® A 


(7) 


C8) 


Tlie pressure ratio corresponding to this I'fach., number ls_ obtained from 
the isentroplc equation ‘ _ 
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E T 


(9) 


"by substituting ©g,uation (4) into equation 


( 9 ) and integrating to give 


( £ e e^/^o 

X = 1 _ e^Ao ee/54-i" Tq ~J/T^ 

Po VTo/ 1 _ 00/^1 ^ 

Similarly, the density ratio is found to be 

(e e^A e e^Ao \ 
p /T\5/2 1 _ e®Ao ©6/^1 To e6/'^o_2_J 
Po “ VToj ;j_ _ 00 /T ® 

The differences inrolTed in the us© of equations (l) and (2) to predict 
the tenrperature, density, and pressure ratios corresponding to a given 
Mach number are given in figures 2(a) and 2(b), in terms of the percentage 
differences from the value given by equations ( 7 ), (lO), and (U) for 
stagnation temperatures of 1000 ° and 2000 ° F absolute. 

It is seen that the assumption of constant heat capacity leads to 
appreciable differences in applying the .isentropic law for a perfect gas 
if sta gn ation temperatures above 1000° F absolute are Involved, 



Errors Involved in the Assumption of the Perfect-Gas 

lav pv = ET for a Gas with Constant Heat Capacities 

In order to evaluate flows in which inqjerfect-gas effects are present, 
an equation of state that takes into account these effects must be chosen. 
For the purposes of this paper an equation which takes into account the 
effects of molectilar forces and size should be sufficient. 

A suitable equation is that of Yan der Waals 


P + 


( 12 ) 
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vhere b is a term correcting for the volume occupied by the molecules 
and a is a term correcting for the effect of molecular forces. 

Figure 3 is a graph of Yan der ¥aals* equation in which the quanti- 
ties p, V, and T have been made nondimionBlonal by dividing by the values 
of these quantities at the critical point Pq^ Vq, and Tg, thus making 
the graph suitable for any gas, (See reference 2.) The graph may be used 
for air if an empirical critical point (pc = 37.2 atm, Tq = 238.5° S' aba., 
Vq = 0.6438 slugs/ft3) is assigned to that mixture of oxygen and nitrogen. 
To give this critical point the values of q., b, and E for air when the 
pressure is measured in pounds per square foot, the specific volume in 
cubic feet per slug and the absolute temperature in degrees Fahrenheit 
are a = 8.78 X 105^ b = 0.654, and E = I716. 

The proper eq^iation for an isentroplc expansion of a real gas is 
(see reference 3) 


dS = Cv dT + 


xoT/YssConstant 


dv 


0 - 


(13) 


which for Yan der Waals* equation becomes 


dE = Cy dT + p dv + -% dv = 0 


(14) 


Equation (l4) may be written as 


dE = Cy dT + d(pv) — V dp + -%■ dv = 0 


and since — v dp = w dw 


dE = Cy dT + d(pv) + ■%• dv + w dw = 0 


(15) 


Assuming constant heat capacity at (jonstant volume and integrating 
equation (I5) giTes __ . . . , . 


E = CyT + V fp — ^ " Constant = E, 


(16) 
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This is then the energy equation for a Yan der Waals gas. ElTiding 
through hy the isentropic speed of sound 


c 


2 _ ^ 
" dp 



(t - h)2 



( 17 ) 


and since 



ET _ £a 
V — h v2, 


then 



2CyTo + ^^o 


ET, 




2CvT - 2v 



•y^BT _ 2st, 
(t - h)^ ^ 



(18) 


The value of v for an isentropic expansion to be placed in 
equation (l8) can be farmed from equation (ih) as follows: 



then 


dT 

T 




T — b • 


( 19 ) 


and if , (ly. is constant 


T = 



( 20 ) 


From equations (l8) and (20), Imowing the stagnation conditions for 
an expansion from Tq to T, the Mach number may be calculated. The 
pressure ratio is then found to be 
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Po 


ET a 


T - 13 ^ 



( 21 ) 


and substituting the value of v from equation (20) we obtain 



Figure h shows the conventional pressure ratio and area 
ratio pw/(pw)j^i plotted against Mach number fcr air starting from 
stagnation conditions of 520° F absolute and various pressures compared 
with the value obtained using constant ratio of heat capacities and the 
perfect-gas law, , _ 

Also shown in figure i)- are the values of pw/(pw) 2^{_2 computed by 
Tsien*s method. It is seen that as the Imperfect— gas effects become 
large it is no longer possible to siDg)lify the analysis by neglecting 

terms containing the squares of — and although Tsien’s results 

pv^ 

are in good agreement at 50 atmospheres when the Van der Naals effect 
is moderate. 

It is interesting to note that the speed of sound in a Tan der Waals 

gas 


c 


2 



-B_\ _ 2a 

Ct/ (t - b)2 V 


( 17 ) 


is not equal to 7ET. The espression for the ratio of specific heats in 
a Van der Waals gas is _ 


7 = 1 + 


E / + a 

^ V pv2 - a + 


V 


( 23 ) 



The value of the heat- capacity ratio 7 is plotted in figure 5 as a function of the preeaure 
for air at room tengierature (520° P ahs.) . 

It must be resnembered that these values of 7 are iinpartant fca* calculating the relationship 
between physical guantities in gas flcfws but that these large changes in 7 have only a small 
effect on the fctrces measured on a body in a wind tTinnel, the Inpcorbant parameter in the case of 
forces being the Mach number in the test section. 


Pcrmulas for the Flow of an Tngerfect Gas with Variable Heat Capacities 
Prcm eq^uatlon (15) substituting the value of Cy from equation (3) gives 


E 


f? + (ef 1 


dT + d(pv) + -^ dv + w dv 

y2 


(24) 


and the energy equation beccraes 


5®^ + — + 2v(^ p — + w^ o Constant 

„ 1 \ v^/ 


(25) 


thus 




-/ 

!= 


2B9 _ 2R6 ^ ^ _ 2a^ 

-\r-i 


Fi 1 ^ 1 

v^ET 

5 , (d'f e® A 

(v - b)® 

2 VJ-V / *„ n2 

(e»A - 1 ) J 



(26) 


VO 
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Tie value (jf V to uo uaed in this equation may oe found from equation fl9) 


To “• t j/rp 
v-^ 


B T 


Sutetltuting the value of ^ from equation (3) Into equation (27) and Integrating yields 

H ' 


■I V — t 5 T , 1 — fl efi f^O a eJ^t^ 

log - — ^ ^ ^ log -pp- + log '— ■ ■ ■ - + J ° sL — ^ 


e^Ao £ e^A* Y 

. , ,)Aoy/^ 1 - e®'" _ \^o e^Ac^r^ ee/T_i/ 

+ c^o - ^ e 


The m’semire ratio cOTreepcnadlng to this s3nrmTisj.Qn is a^in found fre 


BT 

T-l 
* “ “flffl 

i'U 


Vo - t Tq2 


■ I I n I (I 
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The value of the ratio of heat capacities 


7 in this case is 


7 = 



pv^ + a 

5 f±f e^A 

2 \Tj _ ^^2 


(29) 


hut the speed of sound is found from eq.uation ( 17 ) hy substituting the 
value of Cy/R from, equation ( 3 ) to he 



1 

1 . fe'f 

2 W 


v%T _ 2g. 
(t - h)2 V 


Figure 6 shows the conventional pressure ratio and area ratio plotted 
against Mach number for air starting from stagnation canditians of 
2000° F absolute and various pressures compared with the value obtained 
using constant ratio of heat capacities and the perfect— gas lav. 


DISCDBSIOIT 


The foregoing analyses show that the effects of variation of heat 
capacities with ten^ierature do not became important in isentropic expansions 
of air until stagnation ten 3 )eratures of the order of 1000° F absolute are 
encountered. Above 1000° F absolute, however, for accurate analysis this 
variation must be taken into account. In general, it may be stated that 

(e«A _ i)= 

becomes appreciable centered to the number 2.5. 

The effeots_ of Fan der Waals’ forces become Important when either 
the temperature is extremely low for near atmospheric pressures or the 
pressure very high for moderate temperatures. These forces must be 
taken into account when the value of a/v^ becomes appreciable compared 
to the pressure p, or b becomes appreciable compared to v. For air 
these effects are unimportant until stagnation pressures of the order 
of 50 atmospheres at stagnation temperature of 520<> F absolute are 
encountered. 


for diatomic gases these effects are Important when f— 
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Tsien^B method agrees well with the reeultB of this ircvestigation 
up to 50 atmospheres In this caae^ hut it appears that it is not 

poBsihle to neglect the squared terms of ^ and — ^ when the effects 

V T^ip 

of Yan der Waals’ forces hecame appreciable. 


CONCnJSICSS 


Tn many cases found in very high Mach number wird tpnnels and in ^ 
flows of high stagnaticai temperature or pressure, . Imperfect— gas effects 
wTirl the effects of variation of heat capacities may he present. 


For diatomic gases the effect of variation of heat capacities hecomes 


important when^ (— 




- 1 )‘ 


hecomes appreciable compared to 5 / 2 . 


For air these effects become appreciable when stagnation conditions of 
1000 ® F absolute or larger are encountered. 


Imperfect— gas effects become important to gas dynamics when a/v^ 
becomes appreciable compared to the pressure , p or b becomes appreciable 
compared to v. When air is expanded from a stagnation temperature 
of 520® F absolute these effects become important if the stagnation presshres 
are of the order of 50 atmospheres or greater. 

Formulas are presented for handling isentroplc expansions tahtog into 
account these phenomena both separately and at the same time. Tsien’s 
method is found to be applicable for small departures from a perfect gas 
but is not accurate when the effects of Yan der ¥aals* forces become 
appreciable. 


Langley Aeronautical Laboratory 

National Advisory Committee for Aeronautics 
Langley Field, Ya. 
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APEEKDIZ • 

DERIVATIOW (F THE TIERATiaHAL HEAT CAPACITY (F A DIATOMIC GAS 

To ao^rive at the Tlhratlonal heat capacitj of a diatomic gas, the 
individual molecules are treated as linear harmonic oscillators of a 
fundamental frequency and ShrSdinger*s equation is solved for the 
allowahle energy states of such an oscillator. These allowahle states 
are then suhstituted into the equation for the canonical energy distrihu- 
tion and the average energy per particle as a function of the absolute 
temperature is found. This may he differentiated to obtain the 
contribution of the vibrational degrees of freedom of the molecule to 
the heat capacity of the gas at any temperature. 

The average vibrational energy per particle found in this vay is 
(see references 4 and 5) 


1 = ^ + 

2 ehvAr _ 1 

vhere 

h Planck’s constant ■ 

V characteristic frequency of mSlscular vibration 

T absolute tenqierature 

i 

Differentiating to obtain the contribution to the heat capacity of this , 
energy yields 


^ib ^ i M ^hv/kT 

E k BT Vkry ^^hv/kT _ 

hv 

For a particular gas = 0 is a constant and may be determined from 
spectroscopic data. The heat capacity ’at constant pressure is then 

^ = I + = I + (i-f 

E 2 E 2 \tJ ^^e/T _ ^^2 
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The -val-ue of 0 for oxygen is if-OlO.lf and for nitrogen 1b 60l^J^.4 fear 
absolute temperatures measured in degrees Fatoenheit, The value 5526 
may be used for air. 
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Figure 1 .- Tariation of lieat capacity, at c<matant pres sure with ten^ierature as glTen hy eq.uatlQn (4). 
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(■b) To = 2000° F a-bs. 
Figure 2.- Concluded. 
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Ratio of specific heats 
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Pressure p , atm 

Figure 5.- Variatioii of the ratio of epeclflc heats , 7 for a Van der 

Waals gas . == 52Q° F ahs . 
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Figure 6.- PresBure ratio and area ratio plotted agatasi 

Imperfect -gaa effects considered). 


Mach number for air (heat capacltj and 
= 2000° F ats. 







